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Abstract 

Let L be a linear difference operator with polynomial coefficients. 
We consider singularities of L that correspond to roots of the trailing 
(resp. leading) coefhcient of L. We prove that one can effectively 
construct a left multiple with polynomial coefficients L oi L such that 
every singularity of L is a singularity of L that is not apparent. As 
a consequence, if all singularities of L are apparent, then L has a left 
multiple whose trailing and leading coefficients equal 1. 



1 Introduction 



Investigation of singular points (also called singularities) of a linear differ- 
ential or difference operator L gives an opportunity to study singularities of 
solutions of the equation L{y) = without solving this equation. 

* Supported by French-Russian Lyapunov Institute grant 98-03 and by RFBR grant 
04-01-00757. 

t Supported by NSF grant 0098034. 
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In the differential case, take L = Y^i=Q ai{z)d'^ G C[z,9] where d is the 
order of L (so aii[z) ^ 0) and d = d/dz. Let S{L) C C be the set of 
roots of ad{z). The finite singularities of L are the elements of S{L). The 
singularity at oo will not be considered in this paper. If p € S{L) and 
if there exist d linearly independent analytic solutions at z = p then p is 
called an apparent singularity. Suppose L has apparent singularities. The 
question is if it is possible to construct another operator L G C[z,9] of 
higher order such that any solution of L{y) = is a solution of L{y) = 0, 
and S{L) = {p € S{L) | p not apparent}. In the differential case the answer 
is affirmative, see [10]. In this paper we give the affirmative answer to the 
corresponding question for the difference case. 

In the remainder of this paper (except the appendix) only the differ- 
ence case will be considered. The shift operator E acts on functions of the 
complex variable z as Ey{z) = y{z + 1). We consider non-commutative 
operator rings C[2;, £■] and C(z)[£^] (the rings of linear difference operators 
with polynomial and, resp., rational function coefficients over C). Let 

L = ad{z)E'^ + --- + ai{z)E + ao{z)e€.[z,E]. (1) 

Assume that the leading coefficient ad{z) and the trailing coefficient a{i{z) 
are both non-zero, and that 00(2;), . . . , ad(^) do not have a non-constant 
common factor. Set ord L = d. 

Definition 1 A root p of ao{z) is called a t-singularity (a trailing singular- 
ity). A root p of ad{z — d) is called an /-singularity (a leading singularity). 

Definition 2 A right-holomorphic (resp. left-holomorphic^ function is a 
meromorphic function on C that is holomorphic on some right (resp. left) 
half plane. In other words, holomorphic when Rez (resp. —Kez) is suf- 
ficiently large. A half-holomorphic function is a function that is right- or 
left-holomorphic. 

Definition 3 A root p of ao{z) (resp. of ad{z — d)) is called an apparent 
t-(resp. I) -singularity if no right-(resp. left) -holomorphic solution has a pole 
at p. An operator L is a t-(resp. l)-desingularization of L if every mero- 
morphic solution of L is a solution of L, and every t-(resp. I) -singularity of 
L is at-( resp. I )-singularity of L that is not apparent. 

We show that both t- and Z-desingularizations exist. We give algorithms 
t-desing and 1-desing for constructing a f-(resp. Z)-desingularization and 
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algorithm desingboth for constructing a desingularization related to both 
trailing and leading coefficients. 

The above definition of a desingularization is not the same as in [1] (see 
the summary of [1] given in Section 4.2). 

Our approach is based on some specific properties of apparent singular- 
ities that are proved in this paper (Propositions 3,5). 

Besides of a theoretical interest, it is useful to have a desingularization 
L of L for solving the continuation problem. Equation L{y) = can be 
used as a tool to define a sequence or a function. If we know the value 
y{z) at every point 2; of a given strip X < Re z < X + S, where S is larger 
or equal to the order of L then we can find the value of y{z) in the strip 
A — 1 < Re z < X, and then in the strip A — 2 < Re 2; < A — 1, and so on. We 
can keep continuing y{z) to the left in this way except when we encounter 
t-singularities. Similarly, we can continue y{z) to the right except at l- 
singularities. Thus, the singularities of L may present obstacles to continuing 
solutions of L{y) =0. If the singularities are apparent, one can overcome 
those obstacles using L instead of L. Another use of desingularization is the 
following: In the process of continuing sequences to the left (resp. to the 
right) one must always divide by the trailing (resp. leading) coefficient. If 
all f-(resp. Z)-singularities are apparent, then one can avoid such divisions 
by computing a desingularization, see Sect. 4.1 for an application. However, 
there is a price to pay, namely that the order increases. 

In the appendix we give algorithm 9-desing for desingularization in the 
differential case. More general results for the differential case can be found 
in [10], where C[z, d] p] C{z)[d]L is computed. We include this appendix for 
completeness and because the proof is short. 

Acknowledgment. The authors thank the anonymous referees for numer- 
ous valuable and helpful comments. 

2 The sets Cg,a(L) and Rq,aiL) 
2.1 The set of singulcirities 

We consider a linear difference operator (1). The set of trailing resp. leading 
singularities is St{L) = {p £ C|ao(p) = 0} resp. Si{L) = {p G C\a(i{p — 
d) = 0}. The set of singularities is S{L) = 5t(L) U Si {L). 

A point p G C is said to be congruent to a i-(resp. /)-singularity of L if 
ao{p + u) = (resp. ad{p — d — u) = 0) for some 1/ G N. Let L be of the 
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form (1) and cxi, . . . , am be all singularities of L. Set 

l{L) = min{oo, Reo"!, . . . ,Keam}, i^{L) = max{— oo,Reo"i, . . . jRea^}. 

(2) 

So L has no singularity in the half-planes Re^; < l{L) and Re^; > k{L). 
Note that any solution F{z) of L{y) = which is defined and holomor- 

phic on a half-plane Re z > k{L) (or, resp., Re z < i{L)), can be continued to 
a mcromorphic solution defined on C whose poles are congruent to i-(resp. 
/)-singularities. 

Starting from this point until the end of Section 3.2 we will consider only 
the t-singularities (the /-singularities can be handled similarly). 

2.2 Systems of linear relations 

Let Q be a non-empty open subset of C which is stable under E, i.e., z € il, 
implies z + 1 e Q,. For example let U be the right half-plane Re^; > k(L). 
Let If : n ^ Che an arbitrary holomorphic function. We associate to (/? a 
new function whose values are formal Taylor power series in e: 

: n C[[s]], 0{z) = ^^^s\ 

Here e is a new variable, rather than a "small number". Of course, when 
£ G C with \e\ small enough the formal series ip{z) converges and its sum is 
equal to ip{z + e). 

If is a polynomial then we can identify (p{z) = ip{z + e). The operator 
L G C[z, E] has polynomial coefficients. We associate to L the operator 

L = ad{z)E'^ + ■■■ + ao{z) = a^iz + e)E'^ + ■■■ + aQ{z + e) 

which acts on functions ^{z) whose values are formal power series in e. 

The operator L acts also on sequences with values in the field C((e)) of 
formal Laurent series. If a finite sequence fq,fq+i,...,fq+d-i € ^((e)) is 
given for some g G C, then, by using the recurrence given by L, one can 
compute series 

fq-l, fq-2, ■ ■ ■ (3) 

An advantage of L in comparison with L is that neither the leading nor the 
trailing coefficient of L vanishes when z is any complex number. However, 
a series fq-m can turn out to be formal Laurent series for some positive 
integer m even when fq,..., fq+d-i formal Taylor series; if ao{zo) = 
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then the series ao(-2o) is not invertible in the ring C[[e]]. The problem of 
defining complex values of solutions at singularities of L is transformed into 
the problem of obtaining formal Taylor series (truncated power series in 
computations). 

Let a be some i-singularity of L. Now choose g G cr + N for which Req > 
k{L). Let be a function whose values are formal power series in e, and 
suppose that !/($) = 0. If the values of $(z) at points q,q + 1, . . . ,q + d — 1 
are formal Taylor series 

Hq) = Fgfi + Fg^ie + Fg,2s'^ + ■•• 

: ; : ; (4) 

^{q + d-l) = Fq+d-1,0 + Fq+d-l,l£ + Fq+d-l,2£^ + • • • 

with Fij G C then using the equality !/($) = we can compute the formal 
Laurent series $(cr), and each coefficient of this series will be a linear form 
in a finite set of i^ij's. This series can contain negative exponents of e. 

We can find conditions on the coefficients -F^j's in (4) that guarantee that 
$(cr) is a Taylor series. Indeed, if we use the generic power series (4) for this 
computation and if we get some terms with negative exponents of e, then 
after equating their coefficients to zero, we get a system of linear relations. 
This system forms a necessary and sufficient condition that initial conditions 
(4) lead to a formal Taylor series $(cr) when $ satisfies L{^) = 0. This gives 
a finite system of linear relations, denoted as Cq^aiL)- Constructing Cq^^{L) 
can be carried out algorithmically, sec also [2] and [5]. The system Cq^aiL) 
is a necessary and sufficient condition that a function ^{z), having values 
(4) at g, . . . , g + d — 1 and satisfying -£($) = has formal Taylor series values 
at z = cr. 

Example 1 For the operator 

{z - l)zE^ - {3z + 7){z - 3)E + {z + 2){z + 1) 

we can take generic formal Taylor series 

$(4) = F40 + F41S + Fi2e^ + • • • 
$(5) = F5o + F^ie + ^52^^ + ••• 

Then the system C4-i(L) is {2OF40 - 39F50 = 0}. This means that i/$(4) 
and $(5) are as above, and if ^ is a solution of L, then $(—1) is a Taylor 
series if and only if 2OF40 — 39F50 = 0. 
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Proposition 1 With notations as above, if Cq^a{L) ^ then Cq^aiL) con- 
tains a relation f{Fqfi, Fq^i^o, . . . , Fq^d-ifl) = 0, where f is a non-zero lin- 
ear form in variables Fq^Q, Fq+i^o, Fg^d-i,o- 

Proof: Write ^{a) = YX^-n P<j,i£' with F^_n not zero. Now A^^ > since 
Cq^a{L) 7^ 0. Furthermore, AT is the highest pole order of any solution of L 
that has formal Taylor series as initial values at the points q,q -\- 1, . . . ,q -\- 

Now -FCT,-iV can be written as a linear combination of the Fq^^j with 
< i < d and < j. Suppose that some Fq^ij with j > appears with 
non-zero coefficient. Now assign the value 1 for this Fq^ij and the vahic 
for all Fqj^ii^ji with {i' 7^ ii,j)- Then we have obtained a function $ with 
initial values in {0, e-'} at q,q + 1, . . . ,q + d — 1 and a pole of order A'' at a. 
Dividing by , we get initial values in {0, 1} C C[[e]] and a pole of order 
N -\- j &t a which is a contradiction since N -\- j > N. □ 

2.3 The set Rq,a{L) 

Let a be some f-singularity of L and g G a + Z. Let ^{z) be a function 
whose values are rational functions of e (we do not expand them into power 
series yet), and suppose that L{^) = 0. If the values of $(2) at points 
q,q -\- 1, . . . ,q -\- d — 1 are known rational functions of e 

+ + (5) 

then we can compute step-by-step the rational functions 

...,$(g-2),$(g-l),...,$(g + d),$(g + d+l),... (6) 

in particular #(cr) using L. Consider the following d-tuples 

(1,0,0,... ,0),(0,1,0,...,0),... (0,0,0,...,!) (7) 

as d sets of initial values (5). For each set of initial values one obtains a 
sequence of the form (6) so we get rational functions $i(cr), i = 1, . . . ,d. 
Let 

i?,,,(L) = {$i(a),...,$d(c7)}; 
so Rq^a{L) consists of d rational functions of e. 

Proposition 2 With notations as above, ifReq > k{L) and Cq^„{L) = 
then no rational function from Rq^a{L) has a pole at e = 0. 
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Proof: Each of the d-tuples (7) is a d-tuple of formal Taylor scries (whose 
non-constant terms are equal to 0). Any rational function from Rq^^iL) can 
be represented by its formal power series. If Cq^cr{L) = 0, then every such 
series has to be a Taylor series. □ 



Remark 1 Let q & a + N as above. Set q' = a + n + 1 where n G Tj is some 
integer for which there exists no integer i > n with ad{cr + i) =0. Suppose 
that q — q' > 0. If ^ is a solution of L, then ^{q), . . . , ^{q + d — 1) can be 
computed from ^{q'), ■ ■ ■ , ^{q' + — 1) without dividing by e. This implies 
that no rational function from Rq'^a-{L) has a pole at e = 0. 



3 Existence of a desingularization; algorithms 

3.1 The sets Cq^a{L) and Rq^a{L) in the case of appeirent sin- 
gularities 

We will use the following known result (see [6, Proposition 4.4 and Theorem 
4.5]). 

Theorem 1 ( Ramis [9]; Barkatou [3]; Immink [6]) The difference equation 
L{y) = admits linearly independent meromorphic solutions Fi, . . . ,F(i that 
are holomorphic in the right half-plane Rez > k{L). Moreover, for some 
sufficiently large integer N , the sequences {Fj{n)}n>N, j = l,---,d, are 
linearly independent, i.e.. 



Fi{n) 
Fd{n) 



Fi{n + d-l) 
Fa{n + d-l) 



7^0 



(8) 



for all n > N. 



As a consequence of Theorem 1 we get that if L is of the form (1), then 
for any complex number q with Re q large enough, there exist meromorphic 
solutions Fi,... ,Fd that are holomorphic in the half plane Re2; > k(-L), 
such that 

Fi(g) ... Fi{q + d-l) 



Fd{q) ... Fa{q + d-l) 



^0. 



(9) 
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Proposition 3 Let a be an apparent t- singularity of L & C[z, then 
Cq^a{L) = for all q e a + N such that Reg > k(L). 

Proof: Suppose that qe(T + N,Req> k{L) and Cq^a{L) ^ 0. Then for 
any q' such that g' — (7 G N we have Cqi^cr{L) 7^ 0, since if $ is a solution of 
L, then $(g')) • • • > ^{q' + d-l) can be computed from . . . , (^{q + d-l) 
without dividing by e. So wc can assume that Reg is larger than any 
preassigned real number. By Definition 3 the functions Fi, . . . ,Ffi mentioned 
in the consequence of Theorem 1 have no pole at a. If Cg^aiL) 7^ then 
by Proposition 1 there exists a non-trivial relation of the form Mo-^g,o + 
uiFq+ifi + • • • + Ud-iFg^fi-1,0 = in Cq^cr{L) and so the columns of (9) are 
linearly dependent which is a contradiction, because inequality (9) has to 
be valid for all g large enough. □ 

Proposition 4 Let a be a t-singularity of L e C[z,E], then Cq^^{L) = 
for all q E a + N such that Re g > k{L) if and only if a is apparent. 

Proof: It is obvious that if Cq^fy{L) = then a is apparent. The converse 
is true by Proposition 3. □ 

As a consequence of Propositions 2, 3 we get 

Proposition 5 Let a be an apparent t-singularity of L £ C[z, E], q G a+N, 
Reg > k{L). Then no rational function from Rq^cr{L) has a pole at e = 0. 

3.2 t- and desinguleirizations 
Algorithm t -de sing. 

Input: L G CfzjE'] with non-zero E^ coefficient. 
Output: A i-desingularization of L. 

1. Let ao,ad G C[z] be the trailing and leading coefficient of L. 

2. Let n G N be the dispersion of ad, ao, which is the largest integer 
such that Od has some root that equals n plus some root of oq. If such 
n G N does not exist then set n equal 0. 

3. Set L2 := ^L G C{z)[E]. 

4. For i from 1 to n, clear the E'^ coefficient of L2 by subtracting ^E^L 
from L2 where c is the E^ coefficient of L2. 
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5. Let 60 € C[z] be the least common multiple of the denominators of 
the coefficients of L2, and set L3 := 60-^2- 

6. Compute s, t € C[z] for which sao + tbo = gcd(ao, bo). 

7. Output: sL + tL^. 

Here is a Maple implementation of t-desing: 

t_desing := proc(L, E, z) 

local aO, ad, n, L2, i, bO, L3; 

aO := coef f (L,E,0) ; # trailing coefficient 

ad := lcoeff(L,E); # leading coefficient 

n := LREtools [dispersion] (ad, aO, z, 'maJcimalO; 

if not type (n, integer) then n := end if; 

L2 := L/aO; # is in C(z) [E] with trailing coefficient 1 

for i from 1 to n do 

# Simplify L2 and clear its E"i coefficient: 

L2 := collect(L2, E, Normalizer) ; 

L2 := L2 - coeff (L2,E,i) * subs(z = z+i, L/aO) * E~i 
end do; 

# Multiply L2 by its denominator to obtain L3 in C[z,E] 
L3 : = primpart (L2 , E) ; 

bO := coeff (L3, E, 0); # Equals the denominator of L2. 
gcdex(aO,bO,z, 's' , 't') ; # gcd(aO, bO) = s*aO + t*bO 
collect (s*L+t*L3, E, factor) # Output: s*L+t*L3 simplified 
end proc: 

Theorem 2 The algorithm t-desing produces a t-desingularization of L. 

Proof: Let a be an apparent i-singularity. If n is a non-negative integer 
then there exists precisely one operator of the form 

L2 = rn+d{z)E^+'' + ■■■ + rn+i{z)E^+' + 1 G C{z)[E] 

that is right-divisible by L (recall that d = ordL). Algorithm t-desing 
computes this operator L2 where n is the dispersion of and qq. If 
u{z) is a solution of L, then it is also a solution of L2 and so u{z) = 
— J2i=i rn+iiz)u{z + n + i). Plugging in 2; = e -|- cr we get 

d 

u{e + f) = - ^ rn+i{£ + a)u{e + a + n + i). 

i=l 
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By taking the initial values (1, 0, 0, ... , 0), (0, 1,0,..., 0), etc., for the u{e + 
a + n + i) we obtain the rational functions in the set Ra+n+i,(T{L) on the 
left-hand side of this equation, and the rational functions —rn+i{£ + cr) on 
the right-hand side. Since a is apparent, the rational functions in the set 
Ro-+n+i,cT{L') have no pole at e = by Proposition 5 when n is sufficiently 
large (the dispersion is large enough by Remark 1). Hence, the — r^+i have 
no pole at a. Thus, a is not a root of the denominators of L2 G C{z)[E]. The 
least common multiple of these denominators equals the trailing coefficient 
60 of L3, so a is not a root of bo. The trailing coefficient of the output 
is sao -|- tbo = gcd(ao,6o)- Since this is a factor of oq, the i-singularities 
of the output form a subset of the i-singularities of L. And since this gcd 
divides 60, it follows that a is not a f-singularity of the output. The same 
argument applies to every apparent i-singularity a, and hence the output is 
a t-desingularization of L. □ 

Therefore the following theorem (the main theorem for f-singularities) is 
proven: 

Theorem 3 Every L G C[2;, i?] is t-desingularizable (in other words, there 
exists a t-desingularization L of L). 

Example 2 For the operator 

L = {2z-l){z-l)E'^ + (5z-l-9z'^ + 2z^^E + z{l + 2z), 

algorithm t-desing returns: 1/3 {-1 + A z) (5 z — 1 — 9 z^ + 2 z^) E^ + 
^2 _ M ^ + 11/3 + ^ ^3 _ 18 ^4 + 8/3 ^5 j ^2 

+ 1/3 (7 + 4z) {5z-l-9z'^ + 2z^)E + l 

The paper [1] gives a proof that the so-called e-criterion is an alternative 
a way to decide if a dcsingularization exists. The e-criterion is based on a 
construction similar to Rq^a{L)- However, in [1] a different definition of a 
desingularization was used; the definition that we use in this paper is stricter 
(as mentioned in Section 1). 

If the dispersion, the number n in the algorithm, is not positive then 
if follows from Theorem 2 that no i— singularitiy is apparent. This fact 
already followed from the approach in [1], see the summary of [1] given in 
Section 4.2. 

Definition 4 We call L a complete t-desingularization of L if L G C[z,E] 
is right- divisible by L and its trailing coefficient is a non-zero constant. If 
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a complete t-desingularization of L exists, then we say that L is completely 
t-desingularizable . 

Proposition 6 A complete t-desingularization of L exists if and only if all 
t- singularities of L are apparent. 

Proof: Suppose a complete t-desingularization L exists. If F{z) is a right- 
holomorphic solution of L, then it is also a riglit-holomorphic solution of 
L. Then F{z) must be holomorphic since the trailing coefficient of L is a 
non-zero constant. Hence all f-singularities of L are apparent. Conversely, 
if all t-singularities of L are apparent then a complete f-desingularization 
exists by Theorem 2. □ 

Algorithm t -de sing removes at least all apparent t-singularities by The- 
orem 2. If a complete t-desingularization exists, in other words, if all t- 
singularities can be removed, then the above proposition shows that algo- 
rithm t-desing will do so. However, this does not imply that t-desing 
always removes as many t-singularities as possible: 

Examples Let L = {z + 2)'^{z-l)^E-{z + l)z{z-2)^ . The t-singularities 
of L are —1,0,2, none of which are apparent. The application o/ t-desing 
to L gives 

-{z + 3){z + 4fE^ + z{z-2f. 

So one t-singularity disappeared even though no t-singularity was apparent. 
Note that the t-singularity can be removed as well: the operator 

A{z + A)^E^ - 3z{z + 3){z + A)E^ + 3(z + 2){z- ifE + 2{z - 2f 

is right- divisible by L. The t-singularity z — 2 can not be removed, because 
if all t-singularities could be removed then all t-singularities would have to 
be apparent, and this is not the case. 

We also implemented an algorithm that removes all singularities that can 
be removed, by reducing this problem to a linear algebra problem over the 
constants. This implementation is available at: 
http : //www . math . f su . edu/~hoei j /papers/desing/ 

and tends to produce nicer desingularizations than t-desing. Wc used 
t-desing in this paper because it is shorter than the linear algebra based 
desingularization algorithm, and because the proof that all apparent singu- 
larities can be removed is easier with t-desing. 
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Example 4 For the operator L 

{z-3){z-2)E + z{z-l) (10) 

we get Ci^i{L) = C4fi{L) = 0. So L must be completely t-desingularizable; 
algorithm t-desing returns: ^ (5 z — 6) (z — 3) (z — 2)^ {z — 1) 
^ (108 + 106 z + 5 + 39 z^) (z - 3) {z-2)E + l. 

Note that nicer desingularizations are possible, for example (E + l)^ of order 
5, or ~{z + l)E'^ + (17z-29)E^ + (17z + 56)E^ + {-z + b)E + l of order 4. 
The operator L from Example 1 is not completely desingularizable because 
C4 -i(L) ^ there. 

For definiteness, we considered the trailing singularities. If one already 
has an implementation of a t-desingularization algorithm, then one can ob- 
tain an /-desingularization algorithm by changing a small number of lines. 
However, one can avoid this duplication of code because one can reduce l- 
desingularization to t-desingularization (and to get the algorithm 1-desing) 
with the following trick: interchange the roles of the leading and trailing 
coefficient by using the automorphism of C[z, E, E^^] given by z 1-^ — z, 
E^ E-^. This trick was used in an implementation [8] of ds, an old naive 
desingularization algorithm [1]. 

3.3 An operator which is a desingularization of L related to 
both trailing and leading coefficients 

So far we considered mainly trailing apparent singularities. Note that The- 
orem 3 is valid, mutatis mutandis, for leading apparent singularities and a 
desingularization related to the leading coefficient. The prefix "It-" indicates 
that we consider both leading and trailing singularities, just like I- and t- 
indicate leading and trailing. The following theorem generalizes Theorem 3. 

Theorem 4 (The main theorem.) Any L G C[z, is It- desingularizable. 

Proof : Let Rt, Ri G C{z)[E], Lt, Li G C[z, E] be such that Lt = RfL (resp. 
Li = RiL) is a t-(resp. Z)-desingularization of L. Consider the operator 

L = Lt + E"'Li, (11) 

where m = max{l, ord Lt — ord Li + 1}. It is clear that L belongs to C[z, E], 
and L is right-divisible by L since L = {Rt + E'^Ri)L. The operator L is 
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simultaneously a t- and an /-desingularization of L because its i-singularities 
are f-singularities of Lt and its /-singularities are Z-singularities of L/. □ 

Analogously to the f-case an operator L is completely Zf-desingularizable 
if and only if any half-holomorphic solution of L{y) = is holomorphic on 
C. In other words, L is completely Zt-desingularizable if and only if all 
singularities (leading and trailing) are apparent. 

We name desingboth the algorithm given above. Thus, the algorithm 
is given in the following manner: 

• Use algorithms 1-desing and t-desing for constructing Li and L^. 

• Return L = Lt + E'^Li, where m = max{l, ordLt — ordL; + 1}. 

Example 5 The operator {z — 2)E — z has complete It— desingularization 
-SE"^ + 3E -1. 

4 An application, previous work, and a conjecture 
4.1 An application of desingulcirization 

As we have mentioned, it is useful to have a desingularization L of L for 
solving the continuation problem. Below we consider another application of 

the desingularization. 

Suppose that the sequence u{0),u{l), . . . satisfies the following relation 

(1 + 16nfu{n + 2) - (224 + 512n)n(n + 1) - (n + 1)(17 + 16nfu{n) = 

which corresponds to L = (1 + IGzfE^ - {224 + 512z)E - (2 + 1) (17+ 16^)2. 
Assume that u{0), u(l) G Z. By substituting z = 0,1,. .. in the relation we 
find: 

u{2) = 289^i(0) + 224u(l), -^(3) = 736n(0) + 578u(l), . . . 

One sees that u{2),u{3) G Z since we assumed that u{0),u{l) G Z. The 
question is now the following: 

Prove that u{n) G Z for every nonnegative integer n. 

Each time we use L to compute the next term u{n + 2) from the two 
previous terms u{n), u{n+ 1) we perform additions, multiplications, and one 
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division, namely by the leading coefficient of L which is (1 + 16n)^. How 
to prove that this division does not cause u{n + 2) be become a fraction? 
This question becomes easy if we find an Z-desingularization w.r.t. the l- 
singularity of L; if we use the algorithm 1-desing (Sect. 3.2), then it 
produces the following operator: 

Li = E^ + (L - -{z + 11)E - 1(143 + 112z){z + 1). 

Now u{n + 3) can be computed from u{n),u{n + l),u{n + 2) using Li. This 
can only introduce powers of 2 in the denominator of u{n + 3) because L; 
has only powers of 2 in the denominator and has leading coefficient 1. So the 
denominators in the sequence must be powers of 2, but must simultaneously 
be odd (and hence equal to 1) because the leading coefficient (1 + 16n)^ of 
L is always odd. Hence u{n) G Z for all nonnegative integers n. 

4.2 Algorithm ds in [1] 

In this section we will review algorithm ds from [1] and compare it with 
algorithm t-desing. Consider the operator 

L = {z-l){z-2){z + l)E^ + {z^ - 3z^ + 3z + 2)E + z'^{z + 2) (12) 

The trailing coefficient has integer roots z = and z = —2. We would like 
to decide if there is a i-desingularization of L, where the t refers to fact that 
we are only considering the trailing coefficient z^(z + 2) of L. We will start 
with the largest integer root first, so wc first only consider the factor z"^. 
The question is if there exists an operator Rt £ C{z)[E] for which RtL is in 
C[z,£^] and has a trailing coefficient that is not divisible by z'^, or even not 
divisible by z. If such Rt G C{z)[E] exists, then there must exist an Rt with 
the same property in 2;~^C[2;~''"][£^]. One sees this by first replacing each 
coefficient of Rt in C(z) by its series expansion in C({z)) = C[[z]][z~^], and 
then throwing away all non-negative powers of z to end up with an element 
of z~-'^C[2;~^]. So, if an operator Rt with the desired property (that RtL is 
in C[2;,£^] and that the factor of the trailing coefficient has disappeared) 
exists in C(2;)[£^], then there exists an operator Rt in 2;~^C[2;~^][£^] with the 
same property. We can write this Rt as 

TV -1 

Rt = J2nE\ n= cjz^ (13) 

j=0 3=-K 
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for some non-negative integers A'^, K and some Cij G C. 

Suppose that there is a Cij ^ for some j < —2. Then take i minimal 
with this property. Now i can not be because then the traihng coefficient 
of RtL can not be in C[z]. The E'^ coefficient of ViE'^L is rj {z + if{z + 2 + i) 
which has a pole of order > 2 at 2; = since {z + i)'^{z + 2 + i) is not divisible 
by z (this is the reason we started with and not with z + 2). Then the E^ 
coefficient of RtL has a pole of order > 2 as well (the pole in the E^ term 
of riE'^L can not cancel against the £"* terms of vqE^L, . . . , rj_i£'*~^L since 
those terms have pole orders < 2 by the minimality of i). This means that 
RtL is not in C[z,E] which is a contradiction. It follows that there can be 
no j < —2 for which Cij ^ for some thus we can take K = 2 without 
loss of generality. 

Now take iV minimal with ^ 0. Then the leading coefficient of RtL is 
rN{z~l + N){z — 2 + N){z + l + N). This must be in C[z], however, only 
has negative powers oiz. It follows that {z—\+N){z—2+N){z+l+N) must 
be divisible by z. Hence N can be no greater than the largest integer root 
of the leading coefficient, which is 2. So we can take N = 2. In the more 
general situation where we want to eliminate several roots of the trailing 
coefficient at the same time one can use the same argument to show that 
one may assume without loss of generality that the order of Rt is bounded 
by the dispersion (the largest root difference in Z) of the leading and trailing 
coefficient of L. 

We now see that if Rt exists, then we may assume it to be of the form 
in equation (13), and from the preceding wc see that we may also assume 
N = 2 and K = 2. This turns the problem into a finite dimensional system 
of linear equations for the Cy. Solving this system decides whether or not 
the factor z^ can be removed, and if so, how to do this. We find the following 
solution 

-1 9 . 5 2 , 1 

Then RtL is in C[2;,i5^] and has trailing coefficient z + 2. Algorithm ds in 
[1] finds this Rt in a slightly different way. Write i?f = E riE^ = z''^ J2 riE' 
where fj = z'^ri € C[z]. We can now work modulo z"^, so we can view fj as 
an element of C[z]/(z^). Now equate {J2riE^)L to zero modulo z^. This 
leads to a system of linear equations for the fj. This system is already in a 
triangular form so it can be solved quickly by applying Gaussian elimination 
over C[z]/{z'^). If we set fo := 1 then one first finds fi = (5/9)z — 2/3. This 
involves a division in C[z]/{z'^), namely one divides by ao(-z + 1) where 
ao{z) = z^{z + 2) is the trailing coefficient of L. After that one computes 
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f2 = — 2:/12, which involves a division by ao{z + 2). 

After eliminating z'^ we can apply a similar process to the factor z + 2 in 
the trailing coefficient. Then one obtains linear equations over C[2;]/(z + 2) 
instead of over C[z]/{z'^) and one can proceed along the same lines, see [1] 
for details. It turns out that z + 2 can be removed as well. 

Algorithm ds tries to remove one root from the trailing coefficient at a 
time. If a root can not be removed, that is, if we encounter a non-apparent 
t-singularity, then algorithm ds stops, and later roots will not be removed 
even if some of them correspond to apparent t-singularities. The algorithm 
t-desing presented in this paper has the advantage that it removes all ap- 
parent t-singularities, even if there are non-apparent t-singularities between 
them. Furthermore, it is shorter than algorithm ds and does not need to 
compute with roots of cq. Algorithm t-desing is also simpler than differen- 
tial desingularization (see Appendix) since it does not need to know which 
singularities are apparent. 

4.3 A conjecture 

Let L G Q[2,ii^] with leading and trailing coefficient ad{z) and ao{d). Let 
A^i be an integer larger than any integer root of a^iz — d)ao{z) and let N2 be 
an integer smaller than any integer root of ad{z — d)ao{z). Then there exist 
d linearly independent sequences ui,. . . ,Ud ■ Ni + N ^ Q and d linearly 
independent sequences vi, . . . ,Vd '■ N2 — N — > Q that satisfy the recurrence 
given by L. Let P denote the set of all prime numbers p that occur as a 
factor in a denominator of a number that appears in a sequence Uj or Vi. So 
p is not in P if and only if every entry of every Ui and Vi does not have p 
in the denominator. We call L smooth if the set P is finite. Note that the 
set P may depend on the choices of the Ui and Vi, but whether P is finite 
or not only depends on L. 

Conjecture: L is smooth if and only if L is completely Zf— desingularizable. 

The conjecture relates analytic properties to number theoretic proper- 
ties, namely it states that sequence solutions (where we consider sequences 
that extend to the right as well as sequences that extend to the left) have 
only finitely many primes in the denominators iff all half-holomorphic solu- 
tions are holomorphic. 

If L has a complete /f-desingularization L then we can extend sequences 
Ui to the right and sequences Vi to the left with L. This will only introduce 
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finitely many primes in denominators (so L is smooth) since L lias constant 
leading and trailing coefficient. This shows that the conjecture is true in 
one direction. 

Appendix: Differential Case 

In this appendix we prove that the statement in Theorem 3 is also true in 
the differential case. The proof is essentially a desingularization algorithm 
(we name it 9-desing). An implementation can be found at: 
http : // www . math . f su . edu/ ~hoe i j / papers/ des ing/ 
However, the results in this appendix are not new; more general results 
were given by Tsai in [10] (a desingularization is one of the elements in the 
output of Tsai's Weyl Closure algorithm). We include this appendix for 
completeness and because the proof is short since the result is less general 
than [10]. 

Let d the differentiation d/dz. If L e C[z, d] and the coefficients of L do 

not have a non-constant common factor, then a singular point or singularity 
of L is a zero of the leading coefficient of L. If L G C(z)[9] and L is monic 

L = a'^ + a„_i(^)5"-^ + --- + ao(z), (14) 

then singularities are poles of a^'s. We will consider differential operators in 
the form (14). Definitions of regular singularity and, resp., irregular singu- 
larity of a given operator L of the form (14) can be found, e.g., in [7]. A 
point that is not a singularity (neither regular nor irregular) is ordinary. 

Definition 5 Let L G C(2;)[0]. A singularity p of L is apparent if there 
exists an open set U with p & U and a basis of solutions of L(y) = that 
are holomorphic on U . 

Definition 6 Let L G C(2)[9] he monic. Then L is called a desingulariza- 
tion of L if it has L as a right-hand factor, and if for every p G C, if p is 
either an ordinary point of L or an apparent singularity of L, then p is an 
ordinary point of L. 

Proposition 7 Let L £ C{z)[d] be monic, have order n, and let p G C 
The following statements are equivalent. 

a) There exist an open set U withp G U and a basis of solutions yi, ■ ■ ■ ,yn 
of L, holomorphic on U, for which yi vanishes at p with order i — 1. 
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b) p is not an irregular singularity of L, the local exponents at p are 
0, 1, . . . , n — 1, and the formal solutions of L at p are in C[[z — p]]. 

c) p is not a pole of any of the coefficients in 'C-{z) of L (i.e., p is an 
ordinary point of L). 

Proof: a) ^ b) is clear; c) =^ a) is Cauchy's theorem; a proof for b) =^ c) 
can be found in [4, Lemma 9.2]. □ 

Proposition 8 Let L G C(^;)[5] he monic, have order n, and let p G C. 
The following statements are equivalent. 

a) p is either an ordinary point or an apparent singularity of L. 

h) L has n linearly independent solutions in C[[2: — p\\ (this is equiva- 
lent to: p is not an irregular singular point of L, the local exponents 
are non-negative integers and the formal solutions at p do not involve 
logarithms, see [4] for more details). 

c) There exists a monic operator L G C(^)[9] that has L as a right-hand 
factor such that p is an ordinary point of L. 

Proof: a) ^ b) is clear; c) ^ a) follows from Cauchy's theorem. For b) =^ 

c), let m be the highest local exponent. Consider the set V of all integers 
< i < m for which i is not an exponent of L at p. Let Li be an operator 
with the following as basis of solutions: L[{z — p)'), i G V. Let L = LiL. 
Then L satisfies part a) of Proposition 7. □ 

Remark 2 In general, solutions in C[[2— p]] need not he convergent, hut h) 
=^ a) of Proposition 8 says that if all formal solutions of L are in C[[z — p]] 
then they are automatically convergent. 

Theorem 5 A desingularization always exists. 

Proof: Let A be the set of all apparent singularities p E C For p E A, let 
m{p) be the highest exponent at p. Let m be the maximum of all m{p), and 
let e(p) C {0, 1, . . . , m} be the set of exponents of L at p. 

Take p € A ioi which m = m{p). Suppose a desingularization L exists. 
The exponents of L at p must be a subset of the exponents of L at p because 
L is a right-hand factor of L. Hence m is an exponent of L at p, but since 
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p is a ordinary point of L it follows that 0, 1, . . . , m are exponents oi L aX p 
as well, so the order of L must be at least m + 1. 

We will now show that a desingularization L of order m + 1 exists. First, 
construct polynomials yi, . . . , ym+i-n G such that for every p G A and 
every i G {0,1, . . . ,m} \ e{p) there is precisely one yj that vanishes at p 
with order i. Let Li be the monic operator whose solutions are spanned by 
L{yi), ... , L{ym+i-n)- Then LiL satisfies condition a) of Proposition 7 at 
every p G A. Thus, p is an ordinary point of LiL for every p ^ A. However, 
L\L need not satisfy the definition of a desingularization because we may 
have created new apparent singularities. Now LiL has order m + 1 and so 
Li has order m + 1 — n. Write Li = d"^'^^~"' + ^aid'^ with G C{z). 
Take hi € C(z) as follows: hi has no poles in C \ ^4, and bi — at must 
vanish at every p € A with order at least M{p) for some integers M(j)). We 
take these integers M{p) high enough to make sure that L2L is an element 
of C(2;)[5] whose coefficients in C{z) have no pole at any p E A, where 
L2 = - ai)d\ Now let L3 = 9™+^-" + J2 hd' and L = L3L. Since the 
hi have no poles in C \ ^, we see that every ordinary point in C of L is an 
ordinary point of L3 and hence an ordinary point of L. The operators LiL 
and L2L have coefficients that do not have poles at any p e A, hence the 
same is true for L = LiL + L2L, and since L is monic it follows that every 
p G ^4 is an ordinary point of L. □ 



Example 6 Let L be the monic operator with ^;cos(2;) and zsm{z) as basis 
of solutions. Then L has one singularity in C, namely at z = 0, which is 

an apparent singularity with local exponents 1 and 2. Thus, to desingularize 
L we must add a solution with exponent 0. Take y\ = z^ = 1 and compute 
L{yi). The explicit form of L is 



2 2 

-d + i + ^ 

z z^ 



so we find L{yi) = 1 + 2/ z'^. Now let Li be the monic operator with 1 + 2/z^ 
as a basis of solutions, then 

A 

Li = d + 



z{z^ + 2) ■ 



Multiplying L on the left by Li adds a solution (namely yi = z^) to L 
with the missing exponent 0. Then L\L satisfies part b) of Proposition 7. 
Hence z = Q is a regular point of LiL. Unfortunately, Li introduces new 
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singularities, namely at z + 2 = 0. We will illustrate how to remedy this 
with a truncated Laurent series at z = 0. We have 

Li = d+ , / = a + 2z-^ - z + - + • • • 
z{z^ + 2) 2 4 

Since the highest power of z in the denominator of L is z^ , we see that if 
we change any z^-term in Li with i >2 then this can not introduce poles at 
z = in the coefficients of L\L. If we remove these terms from, Li we get 
an operator := d + 2z~^ — z for which the coefficients of L^L still have 
no pole at z = 0. Since L and L3 have no other singularities in C apart 
from z = 0, we see that L := L3L has no other singularities in C either, so 
L must he in Cfz, d\ and indeed: 

L = d^ -zd"^ + 3d- z. 

Note that this desingularization process generally makes the singularity 
at z = cx) worse in the sense that the highest slope in the Newton polygon 
at z = 00 is higher for L3 (and hence for L) than it is for L. In Example 
6 this is unavoidable if we want L to have minimal order, although for 
this particular example a desingularization L = {d^ + 1)^ of non-minimal 
order exists whose singularity at z = 00 is not worse than that of L. Our 
implementation always computes a desingularization of minimal order but 
it is not optimal in the sense that it makes no effort to avoid making the 
singularity at z = 00 worse even in examples where this could be done. 

Definition 7 We call L a complete desingularization of L if L G C[z,9] 
has L as right-hand factor and is monic with respect to d. 

Theorem 5 implies that a complete desingularization of L exists if and 
only if all singularities of L in C are apparent singularities. 
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